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^ ' Abstract. We calculate the neutron matter equation of state at finite temperature based on low- 

Q \ momentum nucleon-nucleon and three-nucleon interactions. Our results are compared to the model- 

■ independent virial equation of state and to variational calculations. We provide a simple estimate 
for the theoretical error, important for extrapolations to astrophysical conditions. 
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O . The nuclear equation of state plays a central role in the physics of neutron stars [1] and 

core-collapse supemovae [2, 3]. Renormalization group methods coupled with effective 
field theory offer the possibility of a new and systematic approach to nuclear matter: For 
low-momentum interactions Vio^k [4] with cutoffs around 2fm^\ the strong short-range 
repulsion in conventional nucleon-nucleon (NN) interactions and the tensor force are 
tamed [5, 6]. At sufficient density, Pauli blocking eliminates the shallow bound states, 
and thus the particle-particle channel becomes perturbative in nuclear matter [5]. In 
addition, the corresponding leading-order chiral three-nucleon (3N) interaction becomes 
^ ' perturbative in light nuclei for A < 2fm^^ [7]. Consequently, the Hartree-Fock (HF) 

C ■ approximation is a good starting point, and perturbation theory (in the sense of a loop 

expansion) around the HF energy becomes tractable [5]. The perturbative character is 
due to a combination of Pauli blocking and an appreciable effective range (see also [8]). 
At finite temperature, the loop expansion around the HF free energy can be realized, 
5^ i based on the work of Kohn, Luttinger and Ward [9, 10], by the perturbative expansion 

of the free energy, where the momentum dependence of the self-energy is treated per- 
turbatively. In this work, we include the first-order NN and 3N contributions, as well 
as anomalous and normal second-order diagrams with NN interactions. The pressure, 
entropy and energy are calculated using standard thermodynamic relations. We use the 
cutoff dependence to provide simple error estimates, and find that the cutoff dependence 
is reduced significantly, when second-order contributions are included. 

We start from the perturbative expansion of the grand-canonical potential T, V), 
where /i is the chemical potential, T the temperature and V the volume: 

^1 = ^0 + ^^l,NN + ^^l,3N + ^2,a + ^^2,n + ■ • ■ ■ (1) 



The non-interacting system is given by Q-q, Q.i = fli,NN + ^^i,3N denotes the first-order 
NN and 3N, and a2 a + ^2. n are the second-order anomalous and normal contributions. 




FIGURE 1. The free energy per particle F /N as a function of density p. The left figure gives the first- 
order NN and 3N contributions with a free single-particle spectrum. Second-order anomalous and normal 
NN contributions are included in the right figure. Our results are compared to the virial equation of state 
(virial) [12] and to the variational calculations of Friedman and Pandharipande (FP) [13]. The virial curve 
ends where the fugacity z = e^^l^ is 0.5. 



The free energy F{N,T,V) is obtained by a Legendre transformation of the grand- 
canonical potential with respect to the chemical potential, F{N,T,V) = f2(/i,r, V) + 
llN, with mean particle number A^. Following Kohn and Luttinger [9], we have 



F{N) = Fo(A^) +ai(jUo) +^^2,n(Aio) + 



1 jdQi^nn/d^y 



+ ..., (2) 



where /Jq is the chemical potential of a non-interacting system with the same density 
p = N/V as the interacting system, = -[dQ.o/dii]^^^, and Fo{N) = f2o(/io) + IM)N 
is the free energy of the non-interacting system. The above expansion ensures that 
the r ^ limit is correctly reproduced [9, 10]. The anomalous second-order diagram 
accounts for perturbative corrections to the free single-particle spectrum. 

Our results [11] for the free energy per particle are shown in Fig. 1 for temperatures 
T = 3MeV, 6MeV and lOMeV, where the low-momentum interaction Viow/t is obtained 
from the Argonne vis potential for a cutoff A = 2.1 fm ^ For the 3N contribution at 
the HF level, we find that only the ci and C3 terms of the long-range 2;r-exchange part 
survive (for details on the 3N interaction, see [5, 7]). For the T = 6MeV results, we 
provide error estimates by varying the cutoff over the range A = 1.9fm^^ (lower curve) 
to A = 2.5fm^^ (upper curve). As expected the error grows with increasing density. 
From Fig. 1, we observe that the equation of state becomes significantly less cutoff 
dependent with the inclusion of the second-order NN contributions. 

Li Fig. 1, we also compare our results for the free energy to the model-independent 
virial equation of state [12] and to the variational calculations of Friedman and Pand- 
haripande [13] (FP, based on the Argonne V14 and a 3N potential). We find a very good 
agreement with the virial free energy, and for the densities in Fig. 1 similar results as FP. 




Our results [11] for the energy per particle are presented in Fig. 2. As for the free 
energy, we observe additional binding and a significantly reduced cutoff dependence 
at second order. In contrast to the variational calculation of FP [13], the low-density 
behavior at second order is in good agreement with the virial equation of state [12]. This 
highlights the importance of a correct finite-temperature treatment of second and higher- 
order contributions. This work is part of a program to improve the nuclear equation of 
state input for astrophysics, and to provide error estimates, for example, for the neutron 
star mass and radius predictions. 
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